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Abstract 

Starting from a standard description of an ideal, isentropic fluid, we de- 
rive the effective theory governing a gapless non-relativistic mode — the sound 
mode. The theory, which is dictated by the requirement of Galilei invari- 
ance, entails the entire set of hydrodynamic equations. The gaplessness of 
the sound mode is explained by identifying it as the Goldstone mode associ- 
ated with the spontaneous breakdown of Galilei invariance. Differences with 
a superfluid are pointed out. 
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The hydrodynamics of an ideal, classical fluid was already well understood in the 19th 
century. The case of isentropic flow, for which the entropy per unit mass is constant, is 
particularly simple. The pressure P is here a function of the mass density p only, and the 
flow is automatically a potential flow. A property of such a fluid is that it supports unat- 
tenuated sound waves, i.e., propagating density oscillations. The waves are unattenuated 
because viscosity and thermal conductivity, which usually serve to dissipate the energy of a 
propagating mode, are absent. Sound waves belong to the class of modes having a gapless 
energy spectrum. The purpose of this essay is to argue that their presence in an isentropic 
fluid is an emergent pTopeitj That is, we do not take the existence of this gapless mode 
for granted, but wish to explain it from an underlying principle, namely that of broken 
symmetry. 

To describe the hydrodynamics of an isentropic fluid, we use Eckart's variational principle 
and start with the Lagrange density 

C = y^r^-pe + e[^oP+V■{p^r)], (1) 

where v is the velocity fleld, p the mass density, and e the internal energy per unit mass. 
For isentropic flow e is a function of p alone. The first and second term in (0) represent the 
kinetic and potential energy density, respectively. The variable 6' is a Lagrange multiplier 
introduced to impose the conservation of mass: 

dop + V- (pv) = 0; (2) 

its dimension is [6] = m^s~^. The variation of ([^) with respect to v yields the equation 

V = V^. (3) 

It shows that, indeed, isentropic flow is automatically a potential flow and it also identifles 
the Lagrange multiplier 6 as the velocity potential. (How to incorporate vortices in this 
variational approach will be discussed below.) With (|]), the Lagrange density (|l|) becomes 

C = -p[do9 + l{V9f + e], (4) 
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where we performed integrations by part. A second field equation can be obtained by varying 
the Lagrange density with respect to p. This yields the Bernoulli equation 

doe + \{yef + h = Q, (5) 

with h = d{pe)/dp the specific enthalpy. On taking the gradient of @, and using that for 
isentropic flow one has the thermodynamic relation 

V/i = ^VP, (6) 

with P = p^de/dp the pressure, one obtains Euler's equation 

^ov + iVv^ + -VP = (7) 
P 

governing the flow of the fluid. 

We next wish to investigate the symmetry content of the theory. Classical hydrodynamics 
has the following invariances: 

• Invariance under space and time translations, — + e^, with = (t,x), p = 
0, 1, 2, 3 and a constant vector. 

• Invariance under global translations of the velocity potential, 9 ^ 9 + a, with a a 
constant. 

• Invariance under Galilei boosts, 

t ^ t' = t, X ^ x' = X - ut; 5o ^ <9o = 5o + u ■ V, V ^ V = V, (8) 

with u a constant velocity. 

According to Noether's theorem, symmetries imply conservation laws. More specifically, the 
above invariances lead to the conservation laws: 

• dotoj + ditij = and Sotoo + ditio = 0, where t^i, (/i, z/ = 0, 1, 2, 3) is the energy- 
momentum tensor, with to* = Pi the momentum density and too = the Hamilton 
density; 
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dojo + diji = 0, with jo the mass density and ji the mass current; 

dodoj + diQij = 0, with goj = —Xjjo + tpj and gij = Xiji — tUj the corresponding 
charge densities and currents. Physically, the conservation dG^i/dt = of the charges 
Coi = / d^x Qoi means that the center of mass of the fluid, X = / d^x^p/M^ with 
M = J d^x p the total mass, moves with constant velocity. 



M- , 
dt 



J d^xp. (9) 



Here, the right-hand side denotes the total momentum of the fluid. 

From the Lagrange density (|) we obtain as explicit form for the various charge densities 
and currents 0: 

dC , , 

J. = = (10b) 



dC 



Uj = C6ij - -T^djO = P6ij + pv.vj (10c) 



Pj (= toj) = -QQ^^i^ = (lOd) 
dC 

tio = QQ:^doe = {n + p)v, (loe) 

^oo) = l^^o^ - ^ = + Pe. (lOf) 

Time derivatives do6 have been eliminated through the field equation (^, so that, for ex- 
ample, 

C^ph-pe=(^p^^- ?j (pe) = P. (11) 

A few remarks are in order. First, the Hamilton density Ti. is the sum of the kinetic and 
potential energy density, as required. Second, the equivalence of the mass current j and the 
momentum density p, which is a hallmark of Galilei invariance, is satisfied by the theory. 



Finally, the set of equations (p!OD constitutes all the equations of hydrodynamics. This brings 
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us to the conclusion that the Lagrange density (^) encodes all the relevant equations for the 
description of an isentropic fluid. 

We next turn to the description of sound waves. We will restrict ourselves to waves of 
small amplitude. These generate only small deviations in the density po and pressure Pq of 
the fluid at rest, so that we can expand the Lagrange density in powers of p — po = p, 
with \p\ << Pq: 

C = -{doe + iv^)(po + p) - eopo - hop - ^Kp^ + O(p'). (12) 

The derivatives (') with respect to p are to be evaluated at p = pq. Since for the system at 
rest 6 is constant, it follows from that ho = 0. If we denote the thermodynamic derivative 
dP/dp by c^, which has the dimension of a velocity squared, the coefficient of the quadratic 
term in p can be written as 

h'o = ^P!> = f (13) 
Po Po 

Apart from an irrelevant constant term (— CqPo) the Lagrange density becomes to this order 

C = -{doe + lv')ipo + p)~^p'. (14) 

^Po 

We next eliminate p by substituting 

p = -p^{doe + y), (15) 

Co 

which follows from expanding the field equation (^. Physically, this equation refiects 
Bernoulli's principle: in regions of rapid fiow, the density p = Po + P and therefore the 
pressure is low. After eliminating p, we obtain a Lagrange density governing the velocity 
potential 9: 

Ces = -Po{do9 + iv^) + ^{do9 + iv2)2. (16) 

ZCq 

This is the standard description of a gapless mode in a non-relativistic context [^0. The 
effective theory we derived here is apart from an irrelevant constant identical to a proposal 
by Takahashi P] which was based on symmetry principles. 



The field equation one obtains for tlie velocity potential 9 from ([T6D is non-linear: 

Po{dle + i^ov^) - pclV ■ V + ipo(5ov2 + V ■ Vv^) = 0. (17) 

The information contained in this equation is nothing more than the conservation of mass 
because 9 was initially introduced in (|l|) as a Lagrange multiplier precisely to enforce this 
conservation law. Indeed, remembering that the combination po(f^o^ + denotes Cq times 
p(t,x) = p(t,x) — po, with po the constant density of the fluid at rest, we see that 
reproduces in this approximation. To simplify (|TT|), we replace po in the first and last 
term by the full density p (which is justified to this order) to arrive at the known 0, but 
unfamiliar field equation 

die - clVH = -9ov2 - iv ■ Vv2 (18) 

governing sound waves in the fluid. If we ignore the non-linear terms, it becomes the more 
familiar wave equation 

die - clV^e = 0, (19) 

implying a gapless linear spectrum, and identifying c, which was introduced via the thermo- 
dynamic derivative dP/dp = c^, as the sound velocity. 

The combination dt^e + |(V6')^ appearing in the description of a non-relativistic gapless 
field is dictated by Galilei invariance. To obtain the transformation property of the velocity 
potential e under a Galilei boost (^ we note that since is a velocity field, V6'(t,x) — »• 
Ve'{t',x.') = V0(t,x) — u. This gives as transformation rule for e 

e{t, x) ^ e'{t', x') = e{t, x) - u • x + /(t), (20) 

with f{t) a yet undetermined function of time. To determine f{t) we note that the factor 
—doe in the Lagrange density (0) is the chemical potential per unit mass. Indeed, using the 
standard definition p = dTi/dp, we find 

p = iv^ + /i = -doe, (21) 
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where in the second equahty we used the field equation (|). This identification fixes the 
transformation rule of —doO: 

- doe{t, x) ^ -d'oe'it', x') = -doe{t, x) - u • v(t, x) + iu' (22) 

and in combination with (|20|), yields for f{t) 

ao/(t) = iu2, or fit) = '^nH (23) 

up to an irrelevant constant. It is easily checked that both the combination do9 + |(V^)^ 
appearing in the effective theory (|l^) as well as the field equation (0) are invariant under 
Galilei boosts. So, contrary to what is sometimes stated in the literature , sound waves are 
invariant under Galilei boosts. The wave equation (|1^) is, of course, not invariant because 
essential non-linear terms are omitted. 

From the effective theory ([T6|) one can again calculate the various Noether charge den- 
sities and currents They are, as might be expected, of the same form as the exact 
expressions (|iy), but now with the approximations 

p^Po-^(5o^^ + iv2) (24) 
Co 

as follows from Eq. ([15|) , 

^-fv^ + ^(p-Po)^ (25) 

and 

P ~ -poidoO + y) ^ clip - po). (26) 

This last equation is consistent with the expression one obtains from directly expanding the 
pressure: -P(p) = Pq + pPq since Pq = and Pq = Cq. 

To recapitulate, we have derived an effective theory describing a massless mode starting 
from the Lagrange density (|^) which entails the complete hydrodynamics of an isentropic 
fluid. We now wish to argue that this massless mode is a Goldstone mode associated 
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with a spontaneously broken symmetry. A first indication in favor of such an interpretation 
follows because the Hamilton density Ti displays a property typical for a system with broken 
symmetry, namely that it is not a function of the velocity potential itself, but is a function 
of the gradient of the field. The energy is minimal if 9 is uniform throughout the sample, 
i.e., there is rigidity A second indication comes from the observation that precisely the 
same effective theory as (|16|) arises in the context of (fermionic and bosonic) superfluids P| , 
and there 9 has been identified as a Goldstone mode. 

Usually the broken symmetry can be identified by the general property that the Gold- 
stone mode is translated under the broken symmetry operations. (There may be other effects 
too, but the translation is always there.) Here, this general characteristic does not uniquely 
identify the broken symmetry because 9 is translated under two symmetry operations. Ac- 
cording to the transformation rule (^) with f{t) given in (^), we have that under a Galilei 
boost 

5^9{t, x) := 9'(t, x) - 9(t, x) = -u ■ X + tu ■ W9(t, x), (27) 

where we took the transformation parameter u infinitesimal small so that quadratic and 
higher powers in u may be ignored. The first term at the right-hand side shows that the 
velocity potential is translated under a Galilei boost. The second symmetry under which 9 
is translated is generated by the total mass M = J d^xp. To see this we first compute from 
(D the conjugate momentum 7Tg of 9 

and note that 9 and p are canonically conjugate 0: 

{e(t,x),p(t,x')} = -5(x-x'), (29) 

where { , } denotes the Poisson bracket. We then use a central result of classical field theory 
stating that the charge Q of a continuous symmetry is the generator of the correspond- 
ing transformations of the fields, x(t, x) x'(t, x). More specifically, for an infinitesimal 
transformation (5^x(t, x) = x'(t, x) — x(^)X) one has 
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5^x(t,x) = -a{x(t,x),Q}, (30) 

with a the transformation parameter. Equation (|29|) now immediately shows that under the 
symmetry generated by the total mass, Q is indeed translated 

b^e{t, x) = -a{e{t, x) , M} = a. (31) 

The point is that the generators of Galilei boosts Goj = / d^x{—Xjjo + tpj) also contain the 
mass density p = jo- It is therefore impossible to distinguish a broken Galilei invariance 
from a broken mass symmetry by considering the algebra alone. 

Let us at this point pause and consider the case of a superfluid. It is well established that 
in the normal-to-superfluid phase transition, the mass symmetry is spontaneously broken. It 
is helpful to recall that this comes about because in a superfluid, which is a quantum system, 
many particles Bose-Einstein condense in a single quantum state. This coherence allows us 



to describe the system by a complex field ip{t, x) = y p{t, x)/m exp[im(p{t, x)//;,] normalized 
such that yields the particle number density p/m of the condensate. The field if, which 
will turn out to describe the Goldstone mode of the broken mass symmetry, is a phase field 
and therefore compact. The field ip is usually referred to as the condensate wavefunction, 
even in the modern literature In our view, this is somewhat misleading. As has 



been stressed by Feynman [0, ip is a classical field describing the coherent behavior of 
many condensed particles in the same way as the classical gauge potential of electrodynamics 
describes the behavior of many photons in a single state. For these classical fields there is 
no probability interpretation as is required for wavefunctions ||T^|. 



It has been argued by Feynman |T2| that the ip-field is governed by a non-relativistic 
('?/'*?/^)^-theory defined by the Lagrange density 

= ihrdoi^ - T^IV^P - ^{mri^ - p,f. (32) 
2m 2po 

The potential energy has its minimum along a circle away from the origin at = Po- This 
implies a spontaneous breakdown of the mass symmetry. The Lagrange density (|32D reduces 
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to the one given in (|l^) when the term —h'^{VpY/'Sm?p is ignored. Using the expression for 
the pressure, cf. (0) 



P 



■ (d ^^d \ ^ 



(pe), (33) 



we find that this term gives the contribution —fi'^{SJ'^p)/Am? to the pressure — the so-called 
quantum pressure. The reason for calling it this way is that it is the only place where 
Planck's constant appears in the equations. To the order in which we are working, it is 
consistent to ignore this term. Because the field equation for ip derived from ( |32D has the 
form of a non-linear Schrodinger equation, we will refer to ■?/' as a Schrodinger field. We 
trust however that the reader realizes that it is a classical field unrelated to a Schrodinger 
wavef unction. 

Let us compare the transformation properties of the Schrodinger field with that of the 
velocity potential 9 of an isentropic fiuid. Under a Galilei boost, tlj{t,x.) transforms as 



x) i)'{t\ x') = exp[i(-u ■ X + \\iH)m/n] ij{t, x). (34) 

With rrvp/h denoting the phase of the Schrodinger field, we see that (p transforms in the 
same way as does the velocity potential. 

Using that the canonical conjugate of the ip-field is vr^ = ihtp*, we easily derive the 
Poisson bracket 

{ipit, x), M} = -i{m/h)ip{t, x). (35) 

This shows that the total mass M generates phase transformations on the field: ipiti ^) ~^ 
tlj'{t,x.) = exp{iam/h)tp{t,x.). The phase of the Schrodinger field is consequently trans- 
lated under the symmetry, just like the velocity potential 6. This transformation property 



identifies (p as the Goldstone mode of the broken mass symmetry. The Poisson bracket (|35D 
also implies that ip and p are canonical conjugate cf. 



{^(t,x),p(t,x')} = -5(x-x'). (36) 
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A similar relation holds for superconductors. On quantizing, the Poisson bracket is replaced 
by a commutator. The Heisenberg uncertainty relation that results for the conjugate pair 
has recently been demonstrated experimentally ||14|| . 

As we remarked above, a necessary condition for the breaking of the mass symmetry 
is the presence of a condensate. Such an intrinsic quantum phenomenon, requiring many 
particles in a single state, has no analog in a classical setting. Hence, the mass symmetry 
cannot be broken in classical hydrodynamics. This leaves us with the second possibility, 
namely that of a spontaneously broken Galilei invariance. The breakdown is a result of the 
presence of a finite mass density. This can be inferred P] from considering the transformation 
of the momentum density p(x) under a Galilei boost: 

5^Pi{t, x) = -Uj{pi{t, x), Goj{t)} = -Uip{t, x) + tu ■ Vpi{t, x), (37) 

or with = UjSf 

5fpi{t, x) = -p{t, x.)5ji + tdjPi{t, x). (38) 

If the mass density p is finite, the left-hand side is non-zero, which is a symmetry-breaking 
condition. 

There is an essential difference between the breaking of Galilei invariance and that of mass 
symmetry. Although both symmetries are Abelian, the mass symmetry is a compact U(l) 
symmetry, whereas the Galilei group is non-compact. More specifically, the transformation 
parameter a of the U(l) group has a finite domain (0 < a < 2tt), while the domain of u, the 
transformation parameter of the Galilei group, is infinite. This is closely connected with the 
impossibility to represent the velocity potential of classical hydrodynamics as the phase of a 
complex field. An immediate physical manifestation of the difference is that a system with 
broken U(l) invariance supports topologically stable vortices, whereas a system with broken 
Galilei invariance does not. This is not to say that vortices are absent in the latter case, it 
merely states that their stability is not guaranteed by topological conservation laws. Closely 
connected to this is that the circulation is not quantized in classical hydrodynamics, which is 
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known to exist in superfluids. Yet, the circulation is conserved also in isentropic fluids. This 
is again not for topological, but for dynamical reasons. More speciflcally, the conservation 
is proven by invoking Euler's equation (|^) as was flrst done for an ideal, incompressible fluid 
by Helmholtz [|I^ and generalized to a compressible fluid by Thomson |]16| . 



The easiest way to observe vortices in a classical fluid is to punch a hole in the bottom 
of the vessel containing the fluid. As the fluid pours out a vortex is formed in the remaining 
fluid — a phenomenon daily observed by people unplugging a sinkhole. Often, as happens in, 
for example, superfluid ^He, the core of a vortex is in the normal state so that the symmetry 
is restored there. In the present context this would mean that inside the vortex core, the 
fluid mass density p is zero. This is indeed what is observed: the vortex core consists of air, 
therefore no fluid is present and p = there. 

In the eye of a tropical cyclone — another example of a vortex, nature does its best to 
restore the Galilei symmetry, record low atmospheric pressures being measured there. (A 
complete restoration would imply the absence of air corresponding to zero pressure.) 

It is customary to incorporate vortices in a potential flow via the introduction of so- 
called Clebsch potentials IIT^ . We will not follow this route, but instead use the principle of 



vortex gauge symmetry ||TB| . In this approach, one introduces a so-called vortex gauge fleld 
9^ = , 0^) in the Lagrange density via minimally coupling to the Goldstone fleld: 



d,e^d^e + el, (39) 



with dfj_ = {do, —V) and 



V X 0P = -2a;, (40) 



so that V X V = 2u} yields (twice) the vorticity uj of the vortex. The object 6^ describing 
the vortex is called vortex gauge fleld []T8[ . The combination dn9 + 9^ is invariant under the 



local gauge transformation 

e{t, x) ^ 9it, x) + a(t, x); C ^ C - d.ait, x), (41) 
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with playing the role of a gauge field. The left-hand side of (^0]) may be thought of as 
defining the "magnetic field" associated with the vortex gauge field = V x 0^. 

For illustrative purposes, let us consider an external, static vortex with circulation F 
located along a line L, which may be closed or infinitely long |T^. Then, uj = ^TS{L), 
where d{L) is a delta function on the line L, 

6,{L) = |^%5(x-y). (42) 

This model with a static, external vortex may be thought of as describing the steady fiow in 
the presence of a vortex pinned to a fixed impurity. The field equation for 6 obtained after 
the substitution ( p9D reads 

doi2 + c^V ■ V = dov^ + iv ■ Vv2 - V ■ (43) 

with /i = —{do6 + ) the chemical potential and v = — 0^ the velocity of the fiow in the 
presence of the vortex. The last term gives a coupling of the velocity field to the "electric 
field" associated with 9f,, 

ep = -ve^ - do0^. (44) 

We note that the field equation (|43| ) is invariant under local vortex gauge transformations. 
Ignoring the higher-order terms and choosing the gauge 6^ = 0, we obtain as equation for 
the fiow in the presence of a static vortex: 

V ■ V = 0, or V ■ {Ve - 0^) = 0, (45) 

which is solved by 

^^(x) = -Jd^y ^(x - y)V ■ e^iy). (46) 
Here, G(x) is the Green function of the Laplace operator 

= (47) 

(27r)3 k2 47r|x| ^ ^ 
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Straightforward manipulations then yield the well-known Biot-Savart law for the velocity 
field in the presence of a static vortex p!5|JT8|] 

v^(x) = — f dyx (48) 



An Jl |x — yP ' 

where the integration is along the vortex. This exemplifies the viability of the vortex gauge 
principle as an alternative to describe vortices in a potential flow. 

Let us continue to study the dynamics of vortices — a subject that has recently received 
considerable attention in the literature Because the vortex motion is determined by 



the flow itself, the vortex can no longer be considered as external. We shall see that the 
non-linear part of the field equation (^3]) becomes relevant here. 

In the absence of external forces, the vortex moves with a constant velocity, say. The 
flow in the presence of a moving vortex can be obtained from the static solution ( ^HD by 
replacing the coordinate x with x — v^t. This implies that 

9oVv = -vl- Vv, = -V(vi.-v,). (49) 

To study sound waves in the presence of a moving vortex, we write the velocity field as 
\{t,x) = Vv(x — vit) + V6'(t,x), with 6 describing small variations around the moving 
vortex solution. Equation ( ^9|) then requires that we write for the chemical potential in (^Sf) 



flit, x) = Vi ■ v,(x - Vit) - doHt, x). (50) 



This leads to the linearized field equation in the frame moving with the vortex |19 



d',e{t, x) - c',v'e{t, x) = -v.(x) ■ vdo[2e{t, x) - e{t, o)], (51) 

describing sound waves in the presence of a moving vortex. We have suppressed the prime on 
the coordinate x' = x — v/^t of the moving frame. In deriving ( pT]) we again used the gauge 
9q = 0, and replaced in the last term v/, with V6'(t, 0). This is justified because the vortex 
is carried by the sound wave and in the moving frame, x = denotes the location of the 
vortex. The first term at the right-hand side in (|5l|) stems from the non-linear term Sqv^ in 
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the general field equation (^31). Thus the non-linearity of sound waves becomes detectable. 
Equation ( pTj ) can be used as a basis to study the scattering of phonons of a free moving 
vortex p9|] . 

So far we have contrasted the breakdown of Galilei invariance (caused by a finite density) 
and mass symmetry (caused by a non-zero condensate). A superfiuid, however, has a finite 
density as well as a non-zero condensate. Both symmetries are therefore broken and we 
expect two different Goldstone modes to be present. This is indeed what is observed in 
superfiuid ''He. The system supports besides first sound, which are the usual density waves 
associated with the breakdown of Galilei invariance, also second sound, or entropy waves. 
The latter mode depends crucially on the presence of the condensate and is the Goldstone 
mode associated with the breakdown of mass symmetry. At the transition point, the second 
sound velocity vanishes, whereas the first sound velocity remains finite. This is as expected 
since only the condensate vanishes at the superfiuid phase transition; the total mass density 
remains finite. 
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